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Abstract We study the radiative heat transfer between a spheroidal metallic nanoparticle and a planar 
metallic sample for near- and far-field distances. In particular, we investigate the shape dependence of the 
heat transfer in the near-field regime. In comparison with spherical particles, the heat transfer typically 
varies by factors between 1/2 and 2 when the particle is deformed such that its volume is kept constant. 
These estimates help to quantify the deviation of the actual heat transfer recorded by a near-field scanning 
thermal microscope from the value provided by a dipole model which assumes a perfectly spherical sensor. 



PACS. 44.40.-|-a Thermal radiation - 78.67.Bf Nanocrystals, nanoparticles, and nanoclusters 
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1 Introduction 

The optical properties of metallic nanoparticles depend 
significantly on their shapes, as has been demonstrated, 
e.g.. for elliptical gold particles attached to the apex of 
fiber-based probes for near- field optical microscopy [I]. 
Similarly, the thermal near-field radiation exchanged be- 
tween a nanometer-sized particle at temperature Tp and a 
closely spaced surface at different temperature Tg is influ- 
enced by the particle's shape, although its linear extension 
may be significantly smaller than the dominant thermal 
wavelength. This shape-dependence occurs even when the 
particle's volume, and thus the total amount of radiating 
matter, is kept constant. In this paper we quantify this 
effect for the case of spheroidal metallic nanoparticles. 

Thermally induced near fields have attracted much ex- 
perimental and theoretical attention in last years 0131 . 
Several experiments and models have been designed in or- 
der to measure and describe the radiative heat transfer 
in the near-field regime, i.e., for distances smaller than 
the thermal wavelength. In this regime tunneling of ther- 
mal photons leads to a magnitude of heat transfer which 
can exceed that achieved by black-body far-field radia- 
tion by several orders of magnitude. Possible applications 
of this phenomenon include, in particular, thermophoto- 
voltaic devices [lElinilTllH] . 

The near-field radiative heat transfer between two di- 
electric bodies has been calculated within the framework 
of Rytov's fluctuational electrodynamics [S] for various ge- 
ometries, including two semi-infinite planar bodies and 
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layered structures (see, e.g., [TU|[TT|[T^[T31[ni[T51[TC|P7] ) . 

two spheres (e.g., '18"19"'20[), a sphere above a semi-infinite 
body (e.g. [T3n21,22..23.24.i25] ). and certain other two-di- 
mensional geometries (26| . Early experiments to detect 
the radiative heat transfer between two effectively semi- 
infinite bodies with flat surfaces have been performed by 
Hargrcaves , and by Domoto et al. [55] Also a pioneer- 
ing, but unconclusive experiment [29] by Xu et al. needs 
to be mentioned. An accurate measurement using glass 
plates separated by a micron-sized gap has been reported 
only recently by Hu et al. |30j Moreover, there now ex- 
ist experimental setups measuring the near-field radiative 
heat transfer between a sphere and a fiat sample [3111321 
133] . involving spheres with radii of some 10 pm. 

Relying on the same basic principle, a near-field scan- 
ning thermal microscope (NSThM) has been developed 
by Kittel et al. [M1I351I35] for recording the radiative heat 
transfer at probe-sample distances even down to some 
nanometers. The sensor of this device consists of the tip of 
a scanning tunneling microscope, functionalized to act as 
a thermocouple, so that here the sensor-sample geometry 
differs significiantly from those geometries for which the 
radiative heat transfer has been calculated exactly. The 
foremost part of such a tip typically has a radius of less 
than 50 nm. 

If one regards this active part of the NSThM sensor 
as a sphere, one may employ the familiar dipole model 
(see, e.g., [I31[^[^[^ [^ ) for describing the results ob- 
tained with such an instrument quantitatively. However, 
the actual shape of the sensor is somewhat prolonged, 
resembling more an ellipsoid than a sphere, and varies 
from specimen to specimen, due to the difficult produc- 
tion process. Such variations of the geometry will have 
consequences for the signal measured with the NSThM. 
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This is what motivates our present investigation of the 
shape-dependence of the near-field radiative heat trans- 
fer: In this work, we study the change of the near-field 
radiative heat transfer in response to shape variations of 
ellipsoidal nanoparticles above a flat sample. Besides help- 
ing one to estimate errors implied by the use of the dipole 
model, there may also be other, more general nanotcch- 
nological applications. 

The paper is organized as follows: In Sec. [5] we briefly 
review the dipole model for calculating the near-field heat 
transfer between a dielectric sphere and a planar sample. 
This step mainly serves to collect the required material 
in a form which can be easily generalized. In Sec. [3] we 
extend this dipole model to the geometry of an ellipsoidal 
particle above a flat sample. The near-field heat transfer 
for this geometry is then calculated explicitly in Sec. SI In 
Sec. Owe discuss the shape dependence of the heat transfer 
as the axes of the ellipsoid are varied, and compare our 
results to the heat transfer between a flat sample and a 
sphere. 



2 Dipole model 

The near-field heat transfer between a spherical particle 
with temperature Tp and a sample with a planar surface 
and temperature Ts can be estimated by means of a sim- 
plified dipole model developed previously by several au- 
thors (e.g., (SHIIlinilllllllMj)- To this end, one first 
considers the thermally fluctuating electric and magnetic 
fields Ef and Hf outside the sample, being generated by 
thermally fiuctuating charges in its interior. These fluctu- 
ating fields consist of a propagating and of an evanescent 
part. The heat flux radiated by the sample is given by the 
mean Poynting vector (Ef x Hf ) and yields the well known 
Kirchhoff-Planck radiation law, to which only the propa- 
gating modes contribute, whereas the evanescent modes 
do not flgure here, since they are bound to the surface of 
the sample. 

Now, if an additional spherical particle with radius R 
is placed at a distance a above the sample, as depicted in 
Fig. [TJ the fluctuating flelds Ef and Hf induce an electric 
dipole moment p and a Foucault current jed within the 
particle. As emphasized by Chapuis et al. [25 j this eddy 
current jod is particularly important for metallic particles. 
It causes losses inside the particle which can be described 
by means of an effective magnetic dipole moment rrioff, 
in analogy to the electric dipole moment p. We assume 
a 3> i?, so that higher multipoles may be neglected; in 
principle, this restriction to sufficiently large distances a 
could be removed by including higher multipoles. Within 
the dipole approximation, the energy absorbed by the par- 
ticle can then be written as 



Ef 



(riicff • Hf 



(1) 



the angular brackets denote an ensemble average. For a 
homogenous and isotropic particle the relations between 
the induced dipole moments and the fields arc given by 




Figure 1. Sketch of a spherical nanoparticle with radius R, 
temperature Tp, and permittivity ep placed at a distance a 
above a planar sample with temperature Tg and permittiv- 
ity £S. 



and 



nicff 



Mo Hf , 



(3) 



where ttp and symbolize the electric and magnetic po- 
larizabilities of the particle; Sq and /iq £^re the permittiv- 
ity and the permeability of the vacuum, respectively. In 
general, the polarizabilities have a directional dependence 
and therefore are described by a tensor. For the highly 
symmetric case of a sphere, this tensor reduces to a scalar 
multiple of the unit tensor, so that the polarizabilities arc 
represented by scalar values in this case. It is essential to 
keep in mind that only non-magnetic materials are taken 
into account here, so that the absorption (rhoff • Hf) is 
solely ascribed to the loss due to eddy currents. 

With the help of Eqs. ([2]) and (jSl) one obtains the ex- 
pression 



duj LO 



Im(ap)eo(|Efr). 



-Im(a«)Mo(|Hf| 



(4) 



£0 ftp Ef 



for the energy transferred from the sample to the parti- 
cle [^157] . Here (|Ef |^)^ denotes the frequency-dependent 
correlation function of the electric field, i.e., the expecta- 
tion value of the product of the Fourier transform Ef(a;) 
of the fluctuating electric flcld and its complex conjugate; 
the correlation function (|Hf|^)tj of the magnetic fleld is 
deflncd analogously. 

Reversely, fluctuating charges inside the sphere of tem- 
perature Tp give rise to energy transfer from the sphere 
to the flat surface. The power absorbed within the sam- 
ple, denoted (i-p_s.s)(Tp), can be calculated in the same 
manner as (Ps->p)(2s) above. The resulting expression 
differs from the expression Q only through the temper- 
ature, so that (Pp-^sXTp) = (-Fs->p>(Tp)- As the energy 
(2) flux (Pp_s.s) is directed oppositely to (Ps^p), the resulting 
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overall heat transfer between the particle and the sample 
has the form [24ll37l 



(5) 



(Ps^p) = (Ps^p)(Ts) - (Ps^p)(Tp 



The correlation functions (|Et|^)(^ and (|Hf in Eq. ([4]) 
are known from fluctuational electrodynamics If one 
assumes that the sample occupies the infinite half space 
z < 0, as in Ref. [2], one gets 



2uj 



(2^)2 4fc,o 



(27r)2 27 



Im(r_L) (ej 



l-ki. 



1 e_L 



(6) 



Im(r||)(f 



^11 



for the electric correlation function, where 

huj 



0{u;,T) 



1 



(7) 



is the Bose function with the inverse temperature /3 
\/{k-QT). Moreover, the unit vectors 



e_L := - (-A:y,fcj;,0)' 

II k^K ^ ' 
ef := (fc^7,fcy7,'«^)* 



(8) 



have been introduced, writing k} — k'^^ ky together with 

kzo — \/kQ ~ k'^ and 7 = -^/k^ — fc^ . As usual, kg ~ uj/c, 
where c is the velocity of light in vacuum. In addition, we 
write 



k^_ 



kzo + kz 
gsfczo - kz 
sskzo + kz 



(9) 



for the Fresncl reflection coefficients, with kz 
The first integral in Eq. (jS]), with k < ko, describes the 
propagating part of the fluctuating field, whereas the sec- 
ond integral with k > fcg describes the evanescent part. 
The expression for (-fff.i-fff is obtained from Eq. ^ 
by interchanging r± o r|| and ^ ^ due to a corre- 
sponding symmetry of the electric and magnetic Green's 
functions. 

Knowing the correlation functions above a half space, 
the mean power absorbed by the particle above a planar 



surface can be calculated from Eq. (j?]), provided the polar- 
izabilities ap and ap are given. In the case of a spherical 
particle of radius R these quantities can be derived from 
Mie scattering theory [3H]- Denoting the particle's rela- 
tive permittivity as ep , and introducing the dimensionless 
variables x = koR and y = y/epkoR, one finds [TOIISB] 



= 27rP- 



; (2ep + 1) (sin(y) - y cos(y)) - y'^ sin(y) 



1) (sin(j/) 
-6) 



ycos{y)) + y'^ sm{y) 



2x^ 

2/2 + 3ycot(y)-3) - — 



(10) 



for .T ^ 1, implying that the particle's radius should be 
small compared to the dominant thermal wavelength, R <^ 
Ath- If we demand even \y\ ^ 1, meaning that the radius 
be smaller than the skin depth at thermal frequencies, the 
above expressions reduce to 



AttR 



gEp 



(fc^^)2 _ 

15 



(11) 
(12) 



Evidently, the expression for ap equals the well known 
Clausius-Mossotti formula. 

Before proceeding, we specify the relevant orders of 
magnitude. For bulk metals the relative permittivity is 
well described by the Drude model [39] 



e = 1 - 



(13) 



For gold at room temperature the plasma frequency is 
Wp = 1.4 X 10^^ s^^, while the relaxation rate figures 
as uJt = 3.3 X 10^"^ s~^. Taking the thermal frequency 
Wth ~ lO^'' s^^ one then finds (\/[eifco)~^ ~ 21 nm, set- 
ting a rough upper limit to the radius of a gold sphere for 
which the expressions (fTT|) and ((T2|) may still hold with 
reasonable accuracy. When treating such minuscule parti- 
cles we have to modify the Drude permittivity (|13p in or- 
der to account for surface scattering, since the bulk mean 
free path of the electrons, which is about 42 nm for gold, 
reaches the same order of magnitude as the spheres' di- 
ameters. For spherical particles of radius R the required 
correction is achieved by the replacement |40| 



3 vp 
4"r 



(14) 



where vp is the Fermi velocity. Taking vp — 1.4 x 10® m/s 
for gold, and setting i? = 20 nm, we obtain ujr ~ 5.3 x 



10 



13 „-l 



less than two times the bulk value. On the other 



hand; size quantization effects become essential only for 
radii below 2 nm, and may therefore be neglected here. 



3 Dipole model for ellipsoidally shaped 
particles 

The near-field radiative heat transfer between a flat sam- 
ple with temperature Tg and an ellipsoidally shaped par- 
ticle with temperature Tp. as sketched in Fig. can be 
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expressed through equations similar to Eqs. (U) and ([S]), 
but instead of employing the expressions ap^^ for the po- 
larizabilities of a sphere, the proper polarizability tensors 
for an ellipsoidal particle have to be taken into aceount. 





Figure 2. Sketch of a spheroidally shaped particle with semi- 
axes _Rs and _Rp. The semi-axis with length _Rp lies along the 
spheroid's axis of rotation. 



Denoting the three semi-axes as Rx, Ry, Rz, we re- 



quire a ^ Taa,yi{Rx,Ry,Rz\ = Rmax to substantiate the 
dipole approximation, while keeping i?max smaller than 
the skin depth, thus allowing for approximately homoge- 
neous fields inside the ellipsoids. It is textbook knowledge 
that the induced dipole moment p of an ellipsoidal par- 
ticle is connected with the incident field Ef through the 
relation Hi] 



Pi = eo Ve 



1 



1 + (ep - l)n. 



■Ef, 



(15) 



with Ve = ^RxRyRz specifying the volume of the par- 
ticle, and ep its permittivity. It is assumed here that the 
components of p and Ef are given in the principal-axis 
system of the ellipsoid. The polarizability tensor then is 
diagonal, with diagonal elements given by 



"1. 



ep-1 



1 + (ep - 1)"^ 



(16) 



The quantities rii are the so-called depolarization coeffi- 
cients. They are written as [H] 



rii 



RxRyRz 



ds- 



1 



(,s + Rf)Ris) 



(17) 



These coeffi- 



with R{s) = ^{s + Rl){.s + R^){s + Rl). 
cients depend only on the shape of the ellipsoid, not on 



its volume. They obey the relations 



3 
i=l 



= 1 and 



> 



(18) 



For the special case of a sphere with Rx = Ry ^ Rz ^ R, 
Eq. dnl) yields ni = n2 = 113 = 1/3, so that Eq. (HH) 
correctly reduces to the Clausius-Mossotti formula pT|) . 
In the case of spheroid, i.e., of a rotational ellipsoid with 
two equal semi- axes, Rx = Ry = Rs and Rz = Rp, the 
coefficients take the form HT] 



ni(e) 
^^3(6) 



"2(e) = -(1 -nsie)), 



hi( i±? 1 - 26 



arctan(e)) 



i?a < Rn 



(19) 



i?s > R 



p ' 



with 



1 - 



1 , Rs> R-o 



(20) 



For calculating the effective magnetic polarizability of 
an ellipsoid, we follow a strategy outlined by Tomchuk and 
Grigorchuk [3D], and use the identity 

Wlm(agj/X0(|fff,i|'>- =Rc/ dV (jed • Ef.ed*).; , (21) 

thus shifting the emphasis from the effective magnetic mo- 
ment nicff to a fluctuating eddy field Ef.cd- This requires 
that the material is non-magnetic, so that nieff is entirely 
due to eddy currents jed- The eddy field obeys the equa- 
tions 



V X Ef^cd('^) = iw/ioHf(w) , 
V • Ef,ed(t^) = ; 



(22) 
(23) 



in principle, the second equation constitutes a boundary 
condition for Ef^cd- Because we assume i?max to be smaller 
than the skin depth, Hf can be considered as constant 
within the volume of the particle, so that the eddy field 
Ef_ed depends linearly on the spatial coordinates. The re- 
sulting relation between Ef^od and Hf reads [10] 



-Bf,ed,x(w) = inJ^O-Rx 



i?2 + R2 R2 



Rl 



(24) 



The other components of the electric eddy field are ob- 
tained by cyclic permutation of the indices. We empha- 
size that the underlying assumption of constant fields Hf 
within the particle is not well fulfilled in cases where the 
particle's characteristic linear dimensions are on the order 
of the skin depth, which is about 21 nm for gold at room 
temperature. In such cases, the formalism outlined here 
may still yield the correct orders of magnitude, but not 
exact numbers. 

Next, one has the identity [ID] 



jod(i^) = weoepEf,od(t^) , 



(25) 
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which relates the induced eddy currents to Ef.cd, with £p 
denoting the imaginary part of the particle's permittivity. 
Thus, for a spheroid the absorbed energy (PT|) is found to 
be 

Re / d^r (jed ■ Ef,ed*)i^ = ^A«oepT|-x 
Jvw, c iU 



i-y\ 



(26) 



The effective polarizability is finally read off by comparing 
this expression with the second term in Eq. Again 
the polarizability tensor has non-zero entries only on its 
diagonal, with imaginary parts given by 



^mo^E^x) = Im(aE,a) = 



i?2 



Rl 



(koRsf 
10 



(27) 



For a sphere with Rg = Rp = R, both expressions give the 
correct imaginary part of the Mie formula ((1 



D. 



- 1 




H H 

lm(a£ X ) + \m(a^ ^ ^ 



Rg/ Rp 

Figure 4. Imaginary parts Im(aE,a;) a-nd Im(aE2) of the 
magnetic polarizabilities (|27p at ojth = 10^* for a 

spheroidal gold particle, and their sum, vs. the ratio Rs/Rp. 
The spheroid's volume Ve is kept constant at that of a sphere 
with radius R = 20 nm; surface scattering is taken into ac- 
count. Data are normalized by the imaginary part of the mag- 
netic polarizability Im(ap) of the sphere. 




1 2 3 4 5 



Figure 3. Imaginary parts Im(aE ,j) and Im(Q|; -) of the elec- 
tric polarizabilities (|16p at ujtb = 10^** s^^ for a spheroidal gold 
particle, and their sum, vs. the ratio Rs/Rp. The spheroid's 
volume Ve is kept constant at that of a sphere with radius 
_R = 20 nm; surface scattering is taken into account. Data are 
normalized by the imaginary part of the electric polarizability 
Im(ap) of the sphere. 



of the permittivity. In order to estimate the resulting error, 
we compare in Fig. [5] the electric absorptivities provided 
by the Drude permittivity with the simple surface correc- 
tion ((Ti)) . and without any correction at all. While the ef- 
fect of the correction is clearly visible, the general trends 
are not altered. Hence, a more refined correction taking 
anisotropy into account should not give drastically differ- 
ent results, at least not in the interval 1/5 < Rg/Rp < 5 
considered. 




In Figs. [3] and S] we plot the imaginary parts of the 
electric and of the magnetic polarizability (fT6)) and (|27p 
at the thermal frequency utth — 10^'^ against the ratio 
Rs/Rp, keeping the volume Ve constant at that of a sphere 
with radius i? = 20 nm. The Drude expression (|13p has 
been taken for the permittivity, with the plasma frequency 
ojp = 1.4 X 10^^ s~^ for gold at room temperature, while 
surface scattering has been taken into account through the 
replacement (HH), inserting the geometric mean (R^RpY^^ 
for R. This replacement is not fully correct for nonspheri- 
cal particles [10] : In principle, one then has different scat- 
tering times for different directions, leading to anisotropy 



1e-27 I ' ' ' ' I 

1 2 3 4 5 

Rs/Rp 

Figure 5. Electric absorptivities for spheroidal gold particles 
with a fixed volume corresponding to that of a sphere with 
radius R = 20 nm, both with (corrected) and without (bulk 
data) surface scattering taken into account. 



The curve progressions observed here can be under- 
stood intuitively. For the electric polarizability depicted 
in Fig. [3] the absorptivity Im(Q;g^) decreases monotoni- 
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cally with increasing ratio Rs/ Rp, whereas the absorptiv- 
ity Inr(a| ^) increases. For Rs/Rp ^ 1 the spheroidal par- 
ticle becomes needle-like along the z-direction; hence, the 
polarizability in z-direction is much greater than that in 
X- or y-direetion. On the other hand, when Rs/Rp » 1 the 
spheroid is pancake-like and lies parallel to the x-j/-plane; 
therefore the polarizability in the directions perpendicular 
to the z-direction dominates. 

In the case of the magnetic polarizability shown in 
Fig. |4] one observes that both Im(af ^) and Im(af ^) go 
to zero in the formal limit Rg/Rp — > 0. This is clear be- 
cause these two quantities represent absorptivities caused 
by eddy currents perpendicular to the x- and to the z- 
direction, respectively. Accordingly, for thin needle-like 
spheroids with Rs/Rp <C 1 both absorptivities disappear, 
because no eddy current can arise then. In the opposite 
case, for Rs/Rp ^ 1, eddy currents can be easily induced 
perpendicular to the z-direction in the pancake-like par- 
ticle, whereas perpendicular to the x-direction no eddy 
currents occur, as before. Hence, for Rs/Rp — >■ oo one ex- 
pects Im(ag ^) — > oo and Im(ag ^) — > 0. Since Im(Q!g ^) is 
positive for all values of Rs/Rp and goes to zero in both 
limiting cases there has to be a maximum, which occurs 
near Rs/Rp = 1, as witnessed by Fig. |31 



4 Radiative heat transfer between an 
ellipsoidal nanoparticle and a planar surface 



In accordance with Sec.[5J the radiative heat transfer from 
a planar sample to an ellipsoidal nanoparticle is expressed, 
in analogy to Eq. (HJ, as 



(Ps^E>(Ts)= / diULu Im(ai,Jeo(|£;f..|')^ 

+Im(aE.j)/^o(|-fff,»|^)c^ , 



where summation over repeated indices is implied. It is 
assumed here that the particle is oriented as in Fig. [51 so 
that the surface normal is parallel to a principal axis of 
the ellipsoid. As discussed before, this form (^5)) is quite 
general, requiring only the specification of the polarizabili- 
ties of the nanoparticle, and of the correlation functions of 
the fluctuating fields above the sample's surface. The ex- 
pression (Pe->s) for the energy flux from the nanoparticle 
back to the sample can again be obtained from Eq. (|28p 
by substituting the temperature Tp for Ts. 

When the specific expressions for a spheroid above a 
planar sample are inserted into Eq. (j28p . the propagating 



modes with k < ko yield the contribution 

K<fen 



Im(ai,J[(l-K|2) + -f (l-|r|i|2)] 



+Im(ai,jf2(l-|r|||2) 



(29) 



+Im(a«J[(l-|r||n + -f(l-|r^n] 

+Im(ai,J^(l-|r^n| . 

This reduces correctly to the known result for a sphere by 
setting Rs = Rp = R: 



9 r°° r d^k- 1 

K < k{-] 



Im(ai + aH)[(l-|r||n + (l-Kn] 



(30) 



Observe that sum of the polarizabilities figures here, be- 
cause eo(|Ef P)w = Aio(|Htp)[^ for the propagating modes. 
Accordingly, the relative size of Im(ap) and Im(ap) de- 
termines whether the radiative heat transfer is dominated 
by the electric or by the magnetic part. From Eqs. ([5^ 
and ([5D1) one can calculate the heat radiated by an el- 
lipsoidal or a spherical nanoparticle in the absence of the 
sample [32] by simply setting \r±\ = |r|| | = and inserting 
the temperature Tp instead of Tg. 

The heat transfer mediated by the evanescent modes 
with K > ko from the sample to the spheroidal particle 
takes the form 



cv,S-»E/ — 



r°° r d^K 0-271 

K>fco 



Iin("E.x)[lm(?'-L) + T2lni(r||)] 



(31) 



-Hlm(aE,jT2lm(r||) 

Kg 

+ Im«J [Im(r||) + j^ln-i{ri_)\ 

+ Im(ai.JpIm(r^)| . 

Again this result leads directly to the corresponding ex- 
pression for a sphere [Mj: 



{Pr 



^d.e(.,rs)fco^y^^ 

K>fe0 



i-2 



Im(r||)] (32) 



4-Im(a^)[lm(r||) + ^° Im(r^)] 

fcg 
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Note that in general one has £o(|Efp)i^ 7^ /Lto(|Hfp)^ for 
the evanescent modes. In contrast to the transfer by prop- 
agating modes, this implies that the magnetic contribution 
to the heat transfer, which is proportional to Im(ap), can 
dominate the electric one [23] even if Im(Q;p) > Im(ap). 

In the quasi-static regime, i.e., for distances a even 
smaller than the substrate skin depth, one has k 3> \/jeffco, 
so that Eq. ([ST]) leads to the approximation 



cv,S~iE/ 



X |Im(aE,^ + Q!E.Jlm(?'||) 



(33) 



Furthermore, using the approximate expressions 

Im(rji) 



ImK).|^ 



(34) 
(35) 



for the imaginary parts of the reflection coefficients, the 
integration over the wave number k can be carried out, 
giving 



(P 



1 



x^Im«, + aij^_^^|,^2a)3 
6ia 



(36) 



by the Drude exprcsson with the plasma frequency of gold, 
and with the relaxation rate adapted to surface scatter- 
ing, as before. Since we are interested in the variation of 
the heat transfer caused by alterations of the shape of the 
nanoparticle, we vary the ratio Rg/Rp in such a way that 
the volume Ve of the respective spheroidal nanoparticle 
equals the volume Vp of a sphere with radius R = 20 nm. 
We estimate that the range 1/5 < -Rs/-Rp < 5 is about the 
largest that is still compatible with the dipole approxima- 
tion, and restrict ourselves to this interval. 
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50 nm 
100 nm 
200 nm 
300 nm 
500 nm 
2000 nm 



Rs/Rp 

Figure 6. Electric contribution to the radiative heat trans- 
fer between a spheroidal gold particle with temperature Tp — 
100 K and a planar gold surface with temperature Ts — 300 K, 
normalized to the corresponding value Pp for a sphere. When 
the ratio Rs/Rp is varied, the particle's volume is kept constant 
at that of a sphere with radius _R = 20 nm. Line styles distin- 
guish data for particle-sample distances ranging from 50 nm to 
2000 nm. 



As expected, in the quasi-static regime one has (P^ 



E 

cv.S- 



^ for the magnetic 



for the electric and (PjJ s_j.e) ( 
contribution; these power laws are the same as those for a 
sphere. Since the sums Im(Q!| 2) and Im(ag ^j+aE z) 

appear in the frequency integral, one may expect that in 
this regime (Pev,s->E) behaves like the imaginary part of 
the sum of the polarizabilities at the thermal frequency 
when Rs and Rp are varied. 

It needs to be stressed, however, that this approximate 
result ((36)) has to be taken with some caution, since the 
required short particle-sample distances may conflict with 
those required by the dipole approximation. 



5 Discussion 

Now we study the shape dependence of the radiative heat 
transfer between a spheroidal nanoparticle and a planar 
sample numerically, on the basis of the full expressions 
and (|31l) . To this end, we assume the temperatures Tp = 
100 K for the particle and Ts = 300 K for the sample. The 
permittivity of the sample is described by the bulk Drude 
model (fO)) with parameters for gold; that of the particle 



In Fig. ini we depict the electric contribution to the 
total heat transfer Pg vs. Rs/Rp for several distances a 
between 50 nm and 2000 nm, normalized to the corre- 
sponding energy transfer Pp between the isochoric sphere 
and the sample. Evidently, for small distances the curves 



of P| 



p actually resemble those of Im(a|; ^ 



Fig. [21 as may have been conjectured from the quasi- 
static approximation p6p . notwithstanding its somewhat 
shaky justification. On the other hand, for relatively large 
distances the progression of P^/Pp is similar to that of 
Im(af 2) alone. This behavior stems from the fact that in 
the latter regime, i.e., at distances such that the propagat- 
ing modes dominate the heat transfer, one has (|i?z,pr ^ 



l-^x,pr| 



pr| 

\Ey^pi\'^), as documented in Fig.[7| In contrast. 



in the near-field regime where the evanescent modes domi- 
nate the transfer one finds (|£^z^ovP) = {\Ex^Cv\'^ + \Ey^cv\'^), 

so that Im(a| ^) and Im(a| ^) contribute in equal mea- 
sure there. 

Analogously, we plot in Fig. [5] the normalized mag- 
netic contribution P^ / Pp to the total heat transfer, for 
distances a extending from 50 nm up to 5000 nm. Similar 
to the electric case, in the near-field regime the graphs 



of P^/Pp resemble that of Im(c 



E.x 



z) previously 



8 



Oliver Huth et al.: Shape-dependence of near-field heat transfer 



J/ 



LU 
V 




log(a / m) 

Figure 7. Ratios {\E.^\'^ + \Ey\^) / {\E^\^) (full line) and {\H^\'^ + 
\Hy\^)l{\H^f-) (dashed line) above a semi-infinite planar gold 
sample at Ts = 300 K, as function of the distance a. 



shown in Fig. HI since (|-ffz,cv| ) = (|-H"x,ev| 
for small distances (see Fig. [7]). On the other hand, for 
relatively large distances a the plots of P^/Pp are quite 
similar to that of Im(ag^) alone, since then (|i?2,prp) <C 

HHx^prl"^ + |-ffy,prP)- 



sample. Here wc find that the curves of the ratio P^/P^ 
exhibit a maximum near Rg/Rp = 1, so that the domi- 
nance of the magnetic part is somewhat less pronounced 
for nanoparticlcs with shapes differing from a sphere. For 
markedly needle-like, prolate particles with a very small 
ratio Rs/Rp, the graphs suggest that the electric contri- 
bution could even dominate the magnetic one, because 
the induction of Foucault currents would be suppressed 
by the needle-like shape. We remark that for non-metallic 
nanoparticlcs and non-metallic samples the contribution 
due to induced eddy currents can be neglected anyway, so 
that for such materials it is always the electric contribu- 
tion which dominates the heat transfer. 




I 
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Figure 8. Magnetic contribution to the radiative heat 
transfer between a spheroidal gold particle with temperature 
Tp = 100 K and a planar gold surface with temperature Ts = 
300 K, normalized to the corresponding value Pp for a sphere. 
When the ratio Rs/Rp is varied, the particle's volume is kept 
constant at that of a sphere with radius R — 20 nm. Line 
styles distinguish data for particle-sample distances ranging 
from 50 nm to 5000 nm. 



Figure 9. Logarithm of the ratio Pe/Pb oi the magnetic to 
the electric contribution to the heat transfer vs. Rs/Rp, for the 
situation considered in Figs. [6]and|8l 



We conclude from Figs. [5] and [5] that the near-field 
heat transfer between a metallic nanoparticlc and a sample 
depends to a sizeable extent on the nanoparticlc 's shape 
even when the total radiating volume is kept constant. In 
the example considered, the electric contribution to the 
heat transfer is enhanced for strongly prolate spheroids 
by a factor of about 10 as compared to a perfect sphere, 
and by a factor of about 4 for strongly oblate ones. On 
the other hand, the magnetic contribution is substantially 
reduced in the needle-like case, but exceeds the sphere 
value by a factor of about 2 for pancake-like particles, 
always assuming the orientation specified in Fig. [5J 



6 Conclusions 



The ratio of the magnetic to the electric contribution 
Pe / is drawn in Fig. 13 again for distances from 50 nm 
to 5000 nm. For all distances considered the magnetic con- 
tribution to the heat transfer is much greater than the 
electric one, unless Rs is much smaller than Rp. This fact 
has been discussed in detail by Chapuis et al. [21] for 
a spherical metallic nanoparticlc above a planar metallic 



We have investigated the radiative heat transfer between 
a metallic spheroidal nanoparticlc and a planar metallic 
probe for small and large distances, on the basis of the an- 
alytical expressions ([^^ and pip for the heat transfer me- 
diated by propagating and evanescent modes, respectively. 
By numerical analysis for the example of a spheroidal gold 
particle above a gold surface we have investigated in detail 
the precise shape dependence of the radiative heat transfer 
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in both the near-field and the far-field regime. Further- 
more, we have derived the approximate expression (j36[) 
for the nonretarded evanescent regime, which provides a 
qualitative understanding of the numerical results. 



0. 
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structure, but again the precise shape of the sensor will 
influence the heat current it records. In view of our model 
calculations, we estimate that the values obtained on the 
grounds of the dipole model under the assumption of an 
ideal spherical sensor may deviate, in the appropriate dis- 
tance regime, by up to an order of magnitude from the 
true values, but probably not by more. 
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Figure 10. Heat transfer between a gold spheroid with tem- 
perature Tp — 100 K and a planar gold surface with tem- 
perature Ts = 300 K separated by 100 nm. The volume of 
the spheroid is kept constant at that of a sphere with radius 
R = 20 nm when R^/Rp is varied. The horizontal line marks 
the value of the heat current between the isochoric sphere and 
the sample. 

Figure [TU] shows the absolute total heat transfer be- 
tween a gold spheroid with temperature Tp — 100 K and 
a planar gold surface with temperature Tg = 300 K for the 
separation a = 100 nm, again fixing the spheroid's volume 
at that of a sphere with radius i? = 20 nm when chang- 
ing the aspect ratio. For Rs/Rp = 1/5 the heat transfer 
between spheroid and sample is only about half of that 
between sphere and sample, whereas for Rs/Rp ~ 5 the 
spheroid-sample transfer is roughly two times as efficient 
as that occurring in the sphere-sample geometry. The very 
fact that there exists such a marked shape dependence 
might be of interest for nanoscale thermal engineering, in- 
sofar as it appears possible to control the amount of heat 
transported at nanoscale distances by carefully designing 
the shapes of both the emitting and the absorbing pieces. 

Our discussion is subject to several restrictions: There 
is the dipole approximation ^ made right at the outset, 
the assumption of constant fields Hf inside the particle 
entering into Eq. ([M)) . and the simplified correction (jl4[) 
for surface scattering. Taken together, these simplifica- 
tions render our analysis approximate, rather than exact, 
although they should still capture the essential physics. 
Further issues that should be investigated in future works 
concern the possible effects of spatial dispersion , and 
of surface roughness [S] . 

With respect to the problem of quantifying the actual 
heat transfer in a near-field scanning thermal microscope, 
our study helps to pin down the error margin. A typ- 
ical NSThM sensor tip is larger than the nanoparticles 
considered in our examples, and does possess an internal 
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